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Abstract 

Bi-local mean field theory is applied to one dimensional quantum liquid with 
long range interaction, which has exact ground state wave function. We obtain 
a mean field solution and an effective action which expresses a long range dynamics. 
Based on them the ground state energy and correlation functions are computed. 
The ground state energy agrees fairly well with the exact value and exponents have 
weaker coupling constant dependence than that of partly known exact value. 







Recently there are much interests in one dimensional quantum liquid with poten- 
tial. Ground state wave function and thermodynamic quantities were found by Sutherland [|1| 
and Calogero0. Ground state wave function is of Jastrow type and becomes simple form 
in certain values of coupling constant. It is possible to compute exponents of correla- 
tion functions for these values of coupling constant. For other general values it has been 
impossible to compute them. Kawakami and YangQ computed recently the exponents 
based on conformal field theory. They show that the system belongs to Luttinger liquid 

iii- 

In the present paper we apply a bi-local mean field theory. Mean field theory has 
been successfully applied to many places. It is not clear if a mean field theory can be 
applied to strongly correlated systems0[§. We apply such a mean field theory that 
includes correlations and compare results with partly known exact values. We see that a 
fair agreement is obtained in the ground state energy but an agreement is not good for 
the exponents of correlation functions if the linearized action is used. 

Lagrangian which discribes interacting fermion system is 



where Ao and A are external electromagnetic potentials and V{x — y) is the two body po- 
tential energy between electrons. We concentrate on the repulsive long range interaction, 
V{x — y) = g/\x — yp. Green's functions are computed from the partition function as 

Z = J V^lj^ViPe- io" I "^"^^ , (2) 

where f3 = 1/kT and t = it in the path integral quantization method. The fermionic 
fields ip{x) and ijj^{x) are anti-commuting c- numbers. We write the interaction term as 

-lV'^(xMl/) ^^'^~^V ^W(^)(2-p(a:-i/)) (3) 

I — y) 

+-i)\x)^{x) i)\y)^{y)p{x - y), 

where p{x — ?/) is a c-number function and a suitable form is chosen depending upon the 
interaction and dimension of space. In the present paper we use p{x — y) = 1. 

Let us describe the system with one spatial dimension here, but its extension is 
straightforward. We rewrite the partition function^ using bi-local auxiliary field U{xi, yi, xq) 
that is local in time coordinate and bi-local in space coordinate. 
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Ceff = [{U{xi,yi;xo)U{yi,Xi;xo) - U{xi,yi; Xo)'ilj\y)'ilj{x) (5) 
-C/(yi,a:i;xo)V''^(a:)V'(?/)} - {U{xi,Xi]Xo)U{yi,yi;xo) 

-U{xi,xi;xQ)ip\y)ip{y) - U{yi,yi;xo)ip\x)ip{x)}]^V{x - y), 

jy — J J)Uq~ I dxodxidyi^^^^^^^^{U{xi,yi]xo)U(yi,xi;xo)-U{xi,xi;xo)U{yi,yv,xo)} 

The partition function goes back to the Eq.(2) by integrating the bi-local field U {xi,yi] Xq) 
in Eq.(4). We write Z, by integrating the fermionic field first, as 



N 



JvUe-^^ff^^\ (6) 

Q-S^ff{U) _ j J dxodxi[Co+J dyiC^ff] ^ 

(A) Stationary phase approximation — The integration of the bi-local field is defined 
around a minimum of the classical action Seff{U), i.e., a mean field Uo{xi,yi) which 
satisfies a self-consistency condition, 



dSeffiU) 



dU 



= 0, 

U=Uo 



{ilJ^{x)il:{y))^o=yo =Uo{xi,yi), (7) 

U=Uo 

in stationary phase approximation. If UQ{xi,yi) is a function of xi — yi, the system is 
invariant under translation and momentum is a good quantum number. The ground state 
in a mean field approximation is defined as a state vector in which one particle states up 
to Fermi momentum are filled. Hence the functions 

Ao = Ai = 0, 

Uoixuyi)^ (8) 

7rr 

where r — x\ — yi, satisfy the self-consistency condition, p is the density and is connected 
with the Fermi momentum hj pp = np and p is the chemical potential. Using the above 
Uo{xi, yi) of Eq.(8), Z is written as Z — ZqZ 

Zq = e-^«^/(^°), (9) 
, VUV^jj^V^jje-^^''°+''''^''^^ 

NZn 



2 



S{U, 7/^^,1/;) = J dxodxi{Co + J dyxC^ff). 
The mean field part, Zq, corresponds to the mean field Hamiltonian, 

dXi[^\x)^^{x) + i>\x)F{pi)^{x)] + Ebackground: (10) 

= -2 1 dr^iUoiryP'"- - UoiO)}, 

dxidyi {\Uo{xi - - Uo{x,x)Uo{y,y)}. 

The second term of the kinetic energy and the background energy are generated from the 
mean field part and have the following forms, 



Fip) 



gip"^ + P%)/4:, \p\<pf, 

9Pf\p\/'^, \v\>Pf, 



background 



I dx,{^^. (11) 

J TT 



In mean field approximation, in which we take only Zq, the energy density, the fermion 
propagator and the density correlation function become as 

fPF dp p^ gpl 1 3 gpl 

J-pp 27r 2m 6 tt DTrm 6 tt 



{p{xi,xo)p{yi,xo)) 



sm ppr 

7ir 

11 1 cos2ppr 



iyi)\xx,XQ)i){yx,XQ)) = , (13) 



27J-2 27r^ r 
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Obviously the self-consistency condition Eq.(7) is satisfied. 

Z in Eq.(9) is due to the fluctuations around the mean held, and gives a correction to 
the above energy density and propagator. We compute their effects in the following using 
local field expansions of the bi-local field. The bi-local field is expanded by a complete 
set of relative coordinates as 

U{xi,yi]Xo) = Uo{xi -yi)[l + J^lni^i - 2/i)«n( — - — ,xo)], (14) 

n=0 ^ 

J dr^^[\Uo{r)\^qn{r)qm{-r) - |t/o(0)rg„(0)g^(0)] = 
go = iCo, qi = iCir; Cq, Ci real, 
q2n — 2n— th real polynomials. 
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Q2n+i — i x (2n + 1— th real polynomials). 

For a computational convenience, we modify the long range g/r"^ potential to a short 
range potential, ge~^^ jr^ and let e — > at the end. We see that some of the local fields 
decouple. This complete set is constructed by Schmidt orthogonalization method starting 
from a set of functions {1, r, r^, ■ ■ ■}. With the above normalization of (inif)^ ^^e mass 
terms become diagonal forms of the local fields a„(a;), and the effective action is written 
as. 



dxx{{— + ao(a:))2 + ax{xf + a-2.{x) 
On 



?0 

1 



+ (- + ao(x))V^(x)Fo(p)V(x) + ai{x)i,\x)Fi{v)i>{x) (15) 
+ Z]?n(0)a„(a;)Co(p) g^(0)a^(a;) 

n=0 m=0 

+2{^ E g/(0)a,(x)}V'^(x)V'(x) - i,\x)-^i,{x) + i,\x)^^{x) + ixi,\x)i,{x) 
The coupling strengths of local fields are given by 

V{r) 



Flip) = dr{Uo{r)e'^'qi{r) - t/o(0)gKO)}- ^ 

Co(p) = I dr|^o(0)r^(e^^'- - 1). (16) 

The low energy and long distance physics are determined by Fi{pf)- They behave for 
/ > 2, as 

Fi{pf) = const e3 (/ > 2), (17) 

when the modified form of the potential is used. Obviously Fi[pf) {I > 2) vanish in the 
small e limit, and ai{x) (/ > 2) decouple. 

Hereafter, we take only ao{x) and ai{x) into account and study the system with them. 
The action then is given by 

S{ao, oi, -0, V'^) = / dx[{— + ao{x)f + ai{xf (18) 

+qlao{x)Coip)aoix) + (- + ao{x))ip\x)Fo{p)ip{x) + ax{x)^p\x)Fi{p)^l){x) 

Qo 

+2{^goao(x)}V'^(x)V'(^) - V^(x)^^(a;) + V't(^)|l^(^) + ^,^l,\x)i^{x)\ 
Uq oxq 2m 
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We consider fluctuations of the fermion field of momentum near the Fermi momenta ip^, 
and use the following linearization approximation : 



Fo(p) = Fo{±Pf) + {ptPf)FI,{±Pf), (19) 
ao(a;)Fo(p) = ao(x)Fo(±pF), 

The chemical potential makes single-body energy vanishes at the Fermi momenta ip^r. 

It is convenient to introduce chiral fields ipL and ipR that have a momentum near pp 
or a momentum near —pp- Two chiral fields are further combined into a two component 
Dirac field and the action is written as 

S^C f dx[\{ ^ lijj ,M ^)Y + \^\{^) + go'«o(a;)Co(p)ao(x) (20) 

J vts l + 2Co{p)qo/{UoFo{pF)) vf 

+i>io{po + do)i> + i>'yi{pi + aM: 

where 

Uo 



ao{x) ={Fo{pF) + 2^^}ao{x), 



ai{x) = Fi{pF)ai{x), 

^0 = FoipF), vi = Fi(pf), dx = dxodx^C-\ C = + ^. 

go rn 

The action is reduced, further to 



dx[-^{aQ{xf + r^ax{xf) + —ff,ui^fp.u + higher derivatives] 

= / c^a;[C^{(aoA)2 + r2(9iA)^ + r2(9o0)' + (9i0)2 + 2(l-r^)9oAai0} + ^(a^0)^ (21) 
J Vq ztt 

+ — (5^0)^ + higher derivatives], 

Uu = d^jfl^ - d^a^, ^^ = 81 + Bf, = df,\{x) + e^ud„(f){x), = v^/vl, 
and the inverse of the Dirac operator is written as, 

(2^1 [7/. (P^ + a^)]~^\y) = e*^(*)-^5'^(^)5'o(x - y)e-^^(S)-T«'^(S), (22) 
So{x - y) = (a;|[7MPl~^l?/)' 75 = -«7o7i = ( J ) • 
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Combining Eq.(21) and Eq.(22), we have 



{cosppiy - x)Tr[7oS'o(x - y)] + sinppiy - x)Tr[7iS'o(x - y)]] 
(p(a;)p(2/))xo=yo = const + Tr(7oS'o(a; - yhoSo{y - x)) 



N 



1 cos2»F(a; — v) 

= const + — -r - ^ rf^, 24 

2iT^{x-yy 2'K^{x-yY 

01 = 2 — ^ + --- — tt(^ 0% 



where 



+ 2nC 



CC/vI + 1/2tx' ^ \vI + 2tiC' 

Exponents a and [3f satisfy Luttinger hquid relations but their g dependence are dif- 
ferent from those obtained by Kawakami and Yang based on conformal field theory. 
The reason for this g dependence is due to the coupling dependent kinetic term, F{pi), 
in Eq.(lO) and (20). This term does not exist and I3f goes to infinity as g goes to infinity 
in short range potential models. 

(B) Self- consistent stationary phase approximation — In the previous part, we see 
that the mean field Uo{xi,yi) of Eq.(8) is a solution of classical equation of motion, 
Eq.(7). However, it does not make the total free energy stationary due to the higher 
order correction as is seen in Eq.(23). Especially the long distance behavior of propagator 
is very different from the free propagator. In this section we find a self-consistent mean 
field that includes the long distance fluctuations, and we compute current correlation 
functions under the self-consistent mean field. The mean field now satisfies 

VU{U - {'4)^)u}e-^'ff^^^ = 0, (25) 
{U) = Uo. 

We start from a modified bi-local mean field Uq which incorpolates the radiative cor- 
rection in the long distance region and compute the correction to the propagator. With a 
suitable Uq, the full propagator agrees to the initial Uq and full self-consistency condition 
is satisfied. We find such Uq first and compute other quantities next. 



6 



We start from the following form of the Uo'. 

( suippr f r \1-0F 



U.{r)^\ ^r.' ' (26) 



which has a fractional power in the long distance region and normalized to the density p at 
the origin, and substitute it to Eq.(16). A parameter is unknown in the beginning and 
is determined from the full self-consistency condition. Another parameter ro is connected 
with a short distance dynamics but is regarded as a constant in this paper. 

We repeat the same procedure as before. Namely we first expand the bi-local field 
with a new complete set of functions {qnir)} which satisfies, 

dr'^[\Mr)\^qn{r)qm{-r-) - ^o(0)'gn(0)g„.(0)] = 5,,^ (27) 



as 



°° xi-\-yi 

U{xi,yi;xo) = Uo{xi -yi)[l + ^a„( — - — ;a;o)g„(xi - yi)]. 

n=0 ^ 

The action is expressed with local fields a„ in the same manner as Eq.(15) with couplings 
F;(p) obtained by replacing C/q with C/q. Fi{pp) behaves as, 

Fiipp)^ const e'^, (/>/?f + ^). (28) 

Small e limit of Fi{pf) depends on the exponent (5f- Fi{pf) vanishes in the small e limit 
if / > + |. Hence the number of local fields that couple with the fermionic system in 
this limit depends on the j3p. It is necessary to treat the effective Lagrangian separately 
depending on the magnitude of fip- 

{i) 1 < < \ '■ For this value oi Pf two local fields clq^x) and ai{x) couple with 
fermion and the others decouple in the e — > hmit. Hence the situation is the same as 
Pf — i case. 

(fi) I < /3f < 2 : In this case three local fields do{x), di{x), and 02 (x) couple. The 
interaction parts are given by, 

/ dxilj\x)Fn{p)i'{x) Qyj (x 

), (29) 



n=0 



F„(p) = 2 J dr{Uo{r)Ur)e''^ - t>o(0)g„(0)} 



V{r) 



4 

The couplings F„(p) are approximated with the Fermi momentum value F„(±pi?) as far 
as the long wave length physics is concerned. The couplings F„(±pi?) are even functions 
if n is even and are odd functions if n is odd, 

Fo{-pF)^FoipF), (30) 
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We find two component field representation of the interaction part as, 

ip{x)'yoFo{pF)^p{x)bo{x) + ^{x)'yiFi{pF)^p{x)bi{x), (31) 



Fo{pF) = jFo{pFy + F2{pF: 



2 



, / X Fo{pf)^ , . F2{pf)^ , s 
bo{x) = - ao{x) + a2{x). 

Fo{pf) Fq[pf) 
With the other component bi{x) and 62(2^) defined as 

h ( \ ~ ( \ h ( \ F2{pf). , . , Fq{pf)^ , . 
bi[x) = ai(x), 62 W = -- pTT — r«o(a;) + -~- — -a2{x), 

Fq[pf) Fo{pf) 

the mass term has the equivalent form as before, 6q + 6f + 6| = Oq + + a|. 

(iii) 2 < (Jf '■ More fields couple with fermionic system. In the long distance region 
interaction terms are written into the two component form as. 



iJj{x)'joFo{pF)'ilj{x)bo{x) +ijj{x)'jiFi{pF)ilj{x)bi{x), (32) 



where. 



FoiPF)^J E ^iW = JE Fi 



n=even 



bo{x) = {^F2n{pF)a2n{x)}, 

Fo[Pf) n 

bl{x) = ~l (E F2n+l{pF)a2n+l{x)}. 

Fo{Pf) n 

Thus the fields 60 (x) and bi{x) couple with the fermion in the long distance region with the 
effective coupling strength Fq{pf) and Fi(pf) and the other fields decouple. Consequently, 
the effective Lagrangian has the equivalent form to that of the previous case. 

Now we solve the full self-consistency condition with the effect of the low energy and 
long wave fiuctuations around the mean field included. We have the self-consistency 
conditions, from Eq.(25), 

(bn) = 0, 

PF^output) = (3F{input), (33) 

We find the exponent /Sf numerically as a fimction of the coupling constant g. The 
exponent u = fSp — ^ thus obtained for Pf^o = 1, 0.1 is shown in Fig.l. 

At the end we compute the ground state energy that includes corrections from the 
low energy and long wave length fiuctuations. Since the partition function at the low 
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temperature is expressed as exp(— /3£^o) with the ground energy Eq, the ground state 
energy is found from the partition function. We find the energy in the lowest order from 
the Zq and its correction from the Z. We use the hnerized form of the Lagrangian. This 
Lagrangian is good in pi < pf- The energy density is expressed as 



The expression of AS is vahd in the low energy region. We specify the region of integration 
^ bi| < Pf, bo I < Cpp. The mean field energy density Sq is given in Fig.2, and the 
energy density Sq + AS from these regions are given in Fig.3. In the same figure, the 
exact energy value obtained by Sutherland is given. Wc sec that our vahic is close to the 
exact value. Thus long wave fiuctuations are responsible to the most part of the ground 
state energy. 

As a conclusion, we see that the one dimensional system with long range inter- 
action can be expressed with the self-consistent bi-local mean field theory. The ground 
state energy agrees well with, but the exponents of the correlation functions deviate from 
the exact value if the linearized form is used. 
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Figure Captions 

Fig.l : The exponent z/ = /3f — 1 is given as a function of the coupling constant. The 
dash line shows the value obtained from conformal field theory and dotts show the self- 
consistent mean field values for pfTq — 0.1 and ppfo — 1-0. The real line shows the simple 
mean field value. 

Fig. 2 : The ground state energy density is given. The real line shows the simple mean 
field value, and the dotts show the self-consistent mean field values for pfTq = 0.1 and 
PfTq = 1.0. The dash line is the exactly known value. 

Fig. 3 : The ground state energy density with fluctuations included are given. The real 
line shows the simple mean field value, and the dotts show the self-consistent mean field 
values for PfTq — 0.1 and Pf^o = 1.0. The dash line is the exactly known value. 
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